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Abstract: Fermionic dark matter interacting with the Standard Model sector through a pseudoscalar 
portal could evade the direct detection constraints while preserving a WIMP miracle. We study the 
LHC constraints on the pseudoscalar production in simplified models with the pseudoscalar either 
dominantly coupled to b quarks or t leptons and explore their implications for the GeV excesses in 
gamma ray observations. We also investigate models with new vector-like fermions that could realize 
the simplified models of pseudoscalar portal dark matter. These models yield new decay channels and 
signatures of vector-like fermions, for instance, bbb, brr, and ttt resonances. Some of the signatures 
have already been strongly constrained by the existing LHC searches and the parameter space fitting 
the gamma ray excess is further restricted. On the other hand, the pure r-rich final state is only 
weakly constrained so far due to the small electroweak production rate. 
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1 Introduction 

The nature of dark matter (DM) and its interactions are among the biggest puzzles in astro-particle 
physics. Attempts to understand DM more directly by looking for signals of its interactions with the 
Standard Model (SM) fields are taking place on a variety of experimental and observational frontiers. 
Given the current lack of any non-gravitational evidence for DM, we need to consider all theoretically 
consistent DM possibilities with distinctive experimental signatures. Such an approach can broaden 
the coverage of search programs by probing new, experimentally unexplored regions of the parameter 
space. 

One way to classify DM models is by the portals that connect DM with the SM particles. Among 
the simplest possibilities, the portal could be a pseudoscalar. Pseudoscalar portal dark matter (PPDM) 
has the appealing feature that the interaction it mediates is suppressed by momentum transfer in the t- 
channel process relevant for direct detection experiments, while preserving the weak-scale cross section 
of the s-channel process, which leads to the right amount of thermal relic abundance. 

An experimentally inspired motivation for PPDM is that it offers a simple explanation for hints 
that have been observed in Fermi-LAT gamma ray data. More specifically, data from the Galactic 
center [1-6] and the newly discovered dwarf galaxy Reticulum II [7, 8] can be explained by annihilation 
of DM particles with the mass of a few ten of GeV to about a few hundred of GeV to various SM 
final states [9-66]. The required cross section of DM annihilation is of the same order of magnitude 
as required to yield the thermal relic abundance matching the observations. 

Gamma ray observations are inherently difficult to interpret in the context of DM, mainly due to 
astrophysical processes that not only contaminate any faint signal but also can mimic the expected 
morphological and spectral signatures. For example, even though the gamma ray emission from the 
Galactic center suggests that the morphology is also consistent with the expectation from DM [4], 
it is possible that millisecond pulsars are responsible for the observed excess [67, 68] (for earlier 
discussions, see [69-71]). However, the situation with dwarf galaxies is different. Dwarf galaxies have 
always been considered as the cleanest sources for DM annihilation gamma ray searches. The recent 
gamma ray observations from Reticulum II [7, 8] (but see also [72]) are thus intriguing as they do 
potentially provide for the first time a hint of a DM signature (something of course that needs to be 
confirmed with future explorations). 
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Annihilation of DM particles in SM channels has many implications for cosmology, the most 
interesting one being the effects of energy injection prior and during the epoch of recombination [73- 
78]. The effective ionization energy from DM annihilation can be constrained by CMB measurements, 
the most recent ones being from the Planck Collaboration [79]. Recent studies [77, 78] find that 
light-lepton final states are excluded for cross sections of order < C>(10“^^cm^/s) for annihilations to 
electron-positron pairs, photons, or VV —>■ 4e, where VV is a pair of intermediate vector bosons, each 
of which decays to an electron-positron pair. Similar constraints on the cross section are obtained for 
muon final states. 

Given the challenges in direct, indirect, and cosmological experiments, it will be interesting to 
explore the power of the Large Hadron Collider (LHC) in testing this simple possibility. Indeed the 
pseudoscalars could be produced copiously at the LHC and lead to different final states and signals 
depending on their couplings. If a pseudoscalar (a) is heavy enough, its effect could be encoded by 
dimension-six operators, which have been considered in [80-85] and motivates collider signals with a 
single (&-)jet and missing transverse energy More recently, the particular simplified model in 

the minimal flavor violating framework with pseudoscalar coupling to the SM fermions proportional 
to the SM Yukawa couplings has been extensively studied in [86-89] . We will consider two other types 
of well-motivated simplified models with the pseudoscalar either dominantly coupled to b quarks or 
T leptons among all SM fermions. Since the pseudoscalar couples to both SM particles and DM, it 
could decay either visibly or invisibly, which calls for different search strategies. One needs to take 
into account the corresponding branching fractions in evaluating the collider sensitivity to the full 
parameter space. 

The simple renormalizable pseudoscalar coupling to the SM fermions is not SM gauge 

invariant. It originates from dimension-five operators iaH with H being the Higgs doublet, 
the left-handed fermion weak doublet and the right-handed fermion singlet. This means that the 
pseudoscalar simplified model is generated by integrating out some heavy fields that couple to both SM 
particles and the pseudoscalar. The coupling between the pseudoscalar and the SM fermions is actually 
set by the mass scale of the heavy degrees of freedom. The UV completion (at the LHC energy scale) 
most studied in the literature is the two Higgs doublet model (2HDM) and its variants [21, 23, 55]. 
Yet there exists another simple class of UV completions for the pseudoscalar simplified model, SM 
augmented by vector-like generations of fermions. We will demonstrate that this class of models 
provides new decay channels and signatures of vector-like quarks and leptons, such as triple 5-jet, 
ferr, and ttt resonances. Some of the signatures have already been probed by current multi-jet or 
multi-lepton LHC searches. These constraints on the mass scale of the vector-like generations restricts 
the pseudoscalar coupling to the SM fermions and thus the parameter space that could explain the 
GeV excess. The pure r-rich final states are only weakly constrained at the moment. 

The paper is organized as follows: in Sec. 2, we discuss two simplified models of PPDM in which 
the pseudoscalar couples dominantly either to 6’s or r’s and existing experimental bounds on these 
models. In Sec. 3, we construct the class of vector-like fermion models that generate the simplified 
models and study their new decay channels and LHC signatures. We conclude in Sec. 4. 


2 Simplified Models of PPDM 

The general simplified model of PPDM is 


UppDM = f 


9xX'y^X + '^9ff'l^f I a, 


( 2 . 1 ) 
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Figure 1. Representative diagrams for production of a pseudoscalar in the first simplified model. 

where we assumed DM is a Dirac fermion denoted by x and a is the pseudoscalar. The second 
term sums over all SM fermions (denoted by /), except neutrinos. Notice that the operators, iafj^f 
(or equivalently iaf^ff^ + cc. in terms of Weyl fermions), summed in the second term are not SM 
gauge invariant. They should be understood as dimension-five operators such as iaH^LT'^/M^ and 
iaH^QL b'^n/M^ generated from integrating out some new heavy degrees of freedom beyond the SM 
of mass scale M*. Thus for the simplified model to be valid, gj v/{;\/2M^) with v = 245 GeV 
is restricted to be 1 (otherwise the higher-dimensional operators start playing an important role 
and can’t be ignored). We choose to study the simplified models rather than the effective contact 
interaction operators. The motivation for this choice is because of potential pitfalls of the latter, for 
instance, the breakdown of the effective operator approach at the LHC for a relatively light mediator 
with mass of a few hundred GeV [90-92]. 

There are many different choices of the couplings between a and the SM fermions. The most 
studied one in the literature is gf = cyf, with c a common factor for all SM fermions and yf the 
SM Yukawa coupling. In this case, the pseudoscalar dominantly couples to the top quarks and its 
dominant production channel is the gluon fusion process similar to that of the Higgs boson. The collider 
phenomenology of this case has been studied in [86-89] . The implications for the GeV excesses has been 
further studied in [89], which showed that current searches already constrain much of the parameter 
space that explain the observed excess and that future searches could cover the full parameter space. 
Here, we focus on the following two benchmark simplified models: 

1- 9b » 5/ (/ 7^ V)] 

2 - 9t > 5 / (/ 7 ^ t). 

The choices represent distinctive collider topologies and search strategies. In the first case, the pseu¬ 
doscalar is produced associated with fj-quarks through representative diagrams in Fig. 1. Notice that 
this set of diagrams is beyond the most studied case with s-channel production of the mediator. In the 
second case, the pseudoscalar is produced through electroweak processes, pp 2 t + a with 

a much smaller rate compared to that of a QCD production channel. Both of the simplihed models 
could be easily realized in weak-scale models. For example, they all could arise in SM augmented by 
vector-like generations of fermions with different representations. In addition, the first choice could 
also be realized in (variants of) Type-H 2HDM with tan/3 > 10 [21, 23, 55] and the second choice 
could be realized in (variants of) Type-HI 2HDM. Below we will discuss the current status of each 
case. 

In the first simplified model, the pseudoscalar might also be produced through gluon fusion with a 
heavy-quark loop (e.g., a bottom or top-quark loop). We focus on the associated production assuming 
that the pseudoscalar-top coupling is negligible, but note that this assumption may not hold in every 
UV completion model, e.g., 2HDM with a small tan,5 < 10. In that case, the constraints we derive 
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Figure 2. A 95% CL lower bound on the pseudoscalar-fe coupling as a function of rua for = gt (left) and 
for g^ — igt (right). In both panels, the blue contours give the dark matter annihilation cross section that 
could explain the GeV gamma ray excess in the la range. 


below should be taken as a conservative estimate in those models. The bottom-loop contribution is 
sub-dominant compared to that of the associated production for values of gb that could be probed at 
8 TeV. 

The pseudoscalar can decay to bb and XX with branching fractions determined by and gb'- 


r(Q ^ xx) 9l 

r(a ^ bb) igl ml 


( 2 . 2 ) 


where we neglected the 6-quark mass and the factor of 3 in the denominator is the color factor. There 
are two possible collider searches sensitive to this simplified model. One is the CMS search for the 
pseudoscalar produced in association with 6-quarks and decaying into a 6-quark pair [93], which is 
most constraining when a has a considerable branching fraction to 6 quarks. If a dominantly decays 
into a pair of DM particles, the most sensitive search could be the one with hard 6-jets and large 
missing transverse energy [85, 94]. The CMS Z{bh)H{imj) boost decision tree analysis could also be 
sensitive to this case [95], which we will leave for future explorations. We focus on > 100 GeV. 
The LHC constraints and future reach for region with ma < 100 GeV could be found in [96]. 

Figure 2 shows the most stringent 95% confidence level (GL) bound on the DM coupling to 6- 
quarks as a function of the pseudoscalar mass from either of the two searches for two representative 
cases g^ = gb (left) and g^ = Agb (right) (at a fixed dark matter mass of 48.7 GeV [46]). We also show 
the contours with the DM annihilation cross section that could explain the GeV gamma ray excess 
in the Icr range [46]. We derive the limits in the following way. First we implement the simplified 
models using FeynRules [97], we then generate parton-level events by feeding the model files into 
MadGraph5 [98], and finally we use Pythia8 [99] for parton showering and fragmentation. 

We find that for gb = g^, a decays to 6 quarks with a branching fraction > 75%, and thus the 
direct search for a decaying into 6-quark pairs is more constraining. For g^ = Agb, we find that as 
a decays mostly to DM particles with a branching fraction ~ 85%, the DM search with hard 6-jets 
is more constraining. We only present result up to TOq = 450 GeV because at energies beyond 450 
GeV current searches are only sensitive to gb greater than one, where the simplified model description 
also breaks down. We also take into account the pseudoscalar branching fractions when deriving 
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limits (instead of assuming a 100 % branching fraction of a into one particular channel, as compared to 
previous studies of this simplihed model). In general, for an arbitrary choice of ((76,5 x)j important 
to consider the branching fractions of the pseudoscalar a and study both types of searches in order to 
derive a robust limit. Note that current searches are only sensitive to about > 3 times the coupling 
gt value required for a DM explanation of the GeV gamma ray excess in most of a’s mass range. 

In the second simplihed model, the pseudoscalar is produced in association with two t’s through 
electroweak processes pp —>■ {Z*,'y*) —)■ 2r + a —>■ 4t or 2t + I^t- An interesting question is whether 
the LEP searches were already sensitive to the production of the pseudoscalar through e+e“ —>■ 2T + a 
process if a is light enough. We hnd that the production rate of this channel is too small to have been 
constrained at LEP. For instance, for ma = 50 GeV and = 1, the production cross section is 13 fb 
at ^/s = 200 GeV. Given an integrated luminosity of 233.4 pb“^ collected per LEP experiment in the 
year 2000 [100], this led to only 3 events per experiment produced in that year. However note that 
even the current LHG searches are not sensitive to this electroweak process. For = 50 GeV and 
Qr = 1, the production cross section of a is 16.6 fb. If a dominantly decays to 2r’s, the most sensitive 
search is the CMS multilepton search [101], which only excludes 30 times the cross section for rua = 50 
GeV. Nevertheless, it is perhaps possible to probe such process through the supersymmetry search 
with two hadronically decaying r’s and missing energy if a decays dominantly to DM particles [102]. 
Unfortunately, currently these searches are only sensitive to a production cross section of order a few 
hundred fb and thus this simplified model is still largely unconstrained. 

If, in the second simplified model, the pseudoscalar has a small but non-negligible coupling to 
other SM fermions, e.g., bottoms, this will open up new production channel of the pseudoscalar such 
as the associated production with b quarks. There are constraints from the MSSM Higgs boson search 
in the tt final state [103, 104], which restricts gt to be below 0.03 for = 200 GeV and below 0.1 
for TTia = 700 GeV assuming 100% branching fraction of a —>■ tt. 

3 New Signatures of UV Completions with Vector-Like Fermions 

Concrete models that give rise to the simplified model in Eq. (2.1) could lead to additional interesting 
experimental signatures. One simple UV completion of the PPDM model, Type-H 2HDM and its 
variants, have been extensively studied in [21, 23, 55]. In this section, we focus on another class of 
models that realizes PPDM: the vector-like fermion models. We demonstrate that this model could 
lead to new interesting vector-like fermion signatures that could be worth dedicated searches in the 
LHC Run 2. 

Consider a generation of vector-like quarks, B', B' with charge (3, l)_i /3 and (3, 1 ) 1/3 under SM 
gauge group SU{3) x SU{2)l x U{1)y. We assume that it only couples to the third generation quarks 
of the SM: 

= {Mb&B' + iYaaB% + + cc.) + (3.1) 

where H is the SM higgs doublet and Q 3 is the third generation SM quark doublet. Integrating out 
the B' held leads to a dimension-hve operator 

i^H^Q 3 b%a + cc., (3.2) 

which gives a {Yay^v/y/2MB)abBb'pj. coupling with v = 245 GeV after electroweak symmetry breaking 
(EWSB). 
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Note that the second operator in Eq. (3.1) opens up a new decay channel of the new quarks, 
B' —> ab with a partial width 


r{B' ab) 


Ygms L ml 
487r \ 


(3.3) 


The partial widths of the other standard decay channels B' —>■ VEt, Zb, hb are proportional to y| with 
a ratio r{B' —>• Wt) : r{B' — >• Zb) : r{B' hb) = 2 : 1 : 1, up to some small phase space corrections. 
For Ya > ys, B' —>■ ab could become the dominant decay channel of the new fermions^This leads 
to an interesting new signature of the vector-like B' quarks: 6 6-jets from pair production of B''s 
with the correct combination of jet triplets with the invariant masses near the B' mass and one of the 
three possible jet pairs within each triplet with the invariant mass close to that of the pseudoscalar 
(as shown in Fig. 3, left). In general, the signature of a 6-jet-rich final state (with more than two 6-jets 
and small amount of missing transverse energy) could arise in a variety of new physics scenarios. In 
addition to the vector-like fermion model we discuss here, it could also arise in several supersymmetric 
scenarios, e.g. i?-parity violating and “stealth” supersymmetry [108-110]. Consequently, it could be 
worthwhile to perform a dedicated multi-6-jet search to probe all these different scenarios. 

There are two current searches sensitive to this topology. One is the jet counting analysis in the 
ATLAS multijet search [111]. The analysis requires a large number of hard jets combined with a 
requirement on the number of 6-tagged jets. In particular, the signal region with rijet > 7, > 80 

GeV and n-b-tags = 2 is most sensitive to the signal we are interested in in the context of this work. 
The 95% CL limit on the number of events that new physics contributes to this region is 38. Another 
sensitive search is the CMS search for i?-parity violating gluino decaying to three jets [112]. This 
search uses jet ensemble technique [113] to look for pair-produced three-jet resonances in multijet 
events from pp collisions. The analysis looks for both light- and heavy-flavor jets. The heavy-flavor 
jet search region requires that each event contains at least one 6-tagged jet, which increases the signal 
significance and thus sets stronger constraint. 

To recast the experimental results, we implement the model described by Eq. 3.1 using Feyn- 
Rules, then feed the model files into MadGraphS to generate the parton-level events, and finally use 
Pythia8 [99] to describe parton shower and fragmentation. In this analysis, we use anti-fej' clustering 
algorithm incorporated in Fast Jet [114] to define jets. The cross section of B' pair production was 
taken from [115, 116], which is an NNLO result. We validate our analysis codes using the sample 
of i?-parity violating gluino decays. The constraints on the vector-like quark from the ATLAS jet 
counting analysis is shown in the right panel of Fig. 3. Assuming 100% branching fraction of a —>■ 66, 
the current constraint excludes the new fermion in the mass range 700-820 GeV, which is comparable 
to the results of searches assuming B' decaying to Wt, Zb or hb [117]. Note that the GMS jet res¬ 
onance search requires a fitting to the shape of the jet triplet invariant mass distribution, which we 
approximated with a mass window selection, yielding comparable, yet slightly weaker bounds on the 
B' mass. 

2/3 of order one contributes negatively to the running of the Higgs quartic coupling and could lead to a potential 
vacuum instability problem which appears at around 10 - 100 TeV [105] while either or Ya of order one have Landau 
poles at above a PeV. The model should be taken as a valid effective field theory at energy scales LHC could probe. 

^In addition, Yukawa couplings between the Higgs boson and the new fermions modify various Higgs couplings and 
electroweak precision observables. These constraints are sensitive to the quantum numbers of the new fermions. For 
instance, adding weak singlet vector-like new fermions are strongly constrained by the Zbb coupling measurements yet 
it could be embedded in custodial models such as Refs. [106, 107]. We want to separate the collider studies from the 
indirect probes and focus on the new collider signatures. We refer the readers to [105] and references therein for a 
detailed discussion of indirect constraints on vector-like fermions. 
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To understand the implications of the constraints for a potential DM explanation of the GeV 
gamma ray excess from the Galactic center and Reticulum II, we also plot on right panel of Fig. 3 
bands yielding the cross section of DM annihilating into bb that could explain the excess [46] . We fixed 
to be Qb or so that a —>■ 66 is the dominant decay channel and = v/Mb, which is equivalent 
to setting YaUs = v^. One could see that the constraints on B' forces the pseudoscalar to be lighter 
than 300 GeV in order to explain the GeV gamma ray excess. Then the light a could be probed by a 
direct mediator search as discussed in Sec. 2. 


6 6 



B'^ ba bbb Simplified Model 



ivb [GeV] 


Figure 3. Left: diagram of B' pair production followed by B' —>■ ab,a — >■ bb decay. Right: bound on the 
vector-like B' quarks decaying through B' —>■ ab, a ^ bb with 100% branching fraction based on the results of 
the jet counting analysis in the multijet search in [111] (red solid curve). The green and blue bands give the 
cross section of DM annihilation fitting the GeV gamma ray excess in the bb channel within the la range. We 
fixed to be gi, (green) or gb/3 (blue) so that a ^ bb is the dominant decay channel and gt, = v/Mb, which 
is equivalent to setting VaJ/s = %/2. 

One could also consider the SM augmented with a vector-like generation of leptons. One simplest 
possibility is L',L' with charge (l,l)±i under the SM gauge group SU{3) x SU(2)b x U{1)y- We 
assume that it only couples to the r leptons in the SM with interactions 

Ant = (MlL'L' + iYauL'TB + ysHU^L' + cc.) + iaxi^x, (3-4) 

with £ being the lepton doublet in the SM. Similar to the vector-like B' model, integrating out the 
heavy leptons leads to an effective coupling between a and r leptons. This also leads to a new decay 
channel of L': L' ra ttt. Gonsidering L' pair production, this leads to a 6r final state with 
the two correct r triplet combinations each forming a resonance in the invariant mass spectrum. In 
addition, one of the three r pair combinations in each of the two triplets should also form an invariant 
mass peak at the a mass. Gurrently, the most sensitive experimental probe of this scenario is the 
GMS multilepton search [101]. More specifically, the search regions with one hadronically decaying r 
and a total number of leptons > 4 are the most sensitive ones to the topology of interest. The results 
are shown in Fig. 4. Assuming 100% branching fraction of a —>■ the current constraint only 
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excludes the new leptons below 230 GeV and it is only sensitive to a part of the parameter space that 
could explain the GeV gamma ray excess. 
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Figure 4. Left: diagram of L' pair production followed by L' —>■ ar, a —>■ rf decay. Right: bound on the 
vector-like L' quarks decaying through L' —>■ ar, a ^ rf with 100% branching ratio based on the results of 
the multilepton search in [101] (orange solid curve). The green and blue bands give the cross section of DM 
annihilation fitting the GeV excess in the rr channel within the Icr range. We fixed to be Qt/S (green) 
or 0.1(7r (blue) so that a —>■ rr is the dominant decay channel and Qr = v/Ml, which is equivalent to setting 
VaJ/3 = %/2. 

Gombining the two simplified models above, one could have SM augmented with vector-like B' 
quarks, as well as vector-like leptons L', both of which couple to the third-generation particles. Pseu¬ 
doscalar a couples to both B' and L'. This could lead to a new possibility of B' decay: B' ^ ab ^ brf 
as shown in the left diagram in Fig. 5. In this case, assuming 100% branching fraction, we find a strong 
constraint on this decay topology from the CMS multilepton search [101], which excludes B' below 
800 GeV. This is shown in the right panel of Fig. 5. We also check the multijet search constraint on 
this topology and find a much weaker limit, as demonstrated in the same figure. 

Depending on the representations of the new fermions, there could be more new signatures of 
vector-like fermions beyond the SM. For instance, in a model with a pair of new quarks with charge 
(3, 2 )i/ 6 and (3, 2)_i/g, the new top-quark partner (T') could decay to a + t with a final state of ti+4b 
jets. This decay channel of T' is studied in [118] that shows that current searches still allow for a 
lighter top quark partner with the mass as low as 500 GeV for a certain range of rua. Again in slightly 
more complicated models, it is possible to have T' decaying to a + t followed by a —> which 

leads to a -I- At final state. The strongest constraint in this case comes from the GMS multilepton 
search [101], which excludes T' with the mass up to about 900 GeV, as shown in Fig. 6. Note that 
both Fig. 5 and Fig. 6 could still be compatible with the GeV gamma ray excess, depending on the 
L'’s parameters. 

It is also of interest to embed a simplified phenomenological model of vector-like fermions in a more 
elaborate EWSB model, such as little Higgs models [119, 120]. For example, it has been shown that 











B'^ ba -> brr Simplified Model 




Figure 5. Left: diagram of B' pair production followed by B' — >■ ab, a — >■ rf decay. Right: bound on the 
vector-like B' quarks with the decay topology B' —>■ ab, a —>■ rf based on the multijet search in [111] (red) and 
multilepton search in [101] (orange). 


T'-> ta -> trr Simplified Model 




Figure 6. Left: diagram of T' pair production followed by T' — >■ at,a ^ rf decay. Right: bound on the 
vector-like T' quarks with the decay topology T' at, a ^ rf based on the multilepton search in [101] 
(orange). 


littlest Higgs model could contain a pseudoscalar with the top-quark partner decaying dominantly 
to it and the top-quark in at least part of the parameter space [121]. More recently, more models 
containing the pseudoscalar and new decay channels of top quark partner in the broad context of 
composite EWSB have been proposed in [118, 122]. 
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Finally, we want to comment that all the collider constraints demonstrated in this section are 
based on the assumption that a particular decay chain has a 100% branching fraction. In principle, 
multiple decay topologies exist for a given model. For instance, in the B' model, since the pseudoscalar 
also couples to DM particles, B' could decay through B' ^ ba ^ b + We could then have the 
asymmetric final states with B' —)■ bbb in one decay chain and B' —>■ b + I^t in the other. It is 
rather useful to present these search results in a triangle with Br(i3' —>■ Wt, Zb^ hb), Br{B' —>■ 3b), and 
Bi{B' —)■ b + I^t) =1 at the three vertices. At each vertex, the most sensitive search is quite different, 
as we demonstrate in Fig. 7. 


Br{B' Wt) = 0.5, Br(Z6)=Br(h6)=0.25 



Bt{B' 3b)=l Br{B' 6+MET)=l 

Figure 7 . New triangle of B' search for model described by Eq. 3.1. 

4 Conclusions 

In this paper, we explore two simplified models of pseudoscalar portal dark matter and a class of its UV 
completion with vector-like fermions. We showed that for pseudoscalar dominantly coupled to 6-quarks, 
the most sensitive searches could be either the direct search for scalar particles decaying into 6-quark 
pairs or the 6-jets-|-.^r searches. The current sensitivity still allows for a dark matter explanation of the 
GeV gamma ray excesses in the Galactic center and/or Reticulum II. If the pseudoscalar dominantly 
couples to r’s, current LHC searches are not sensitive to it due to the small electroweak production 
rate. The pseudoscalar coupling to SM could be realized in SM augmented by vector-like fermions 
with the coupling strength set by the mass scale of the new fermions. These new vector-like models 
give rise to interesting signatures such as 666,6rr, rrr resonances. Considering pair production of 
the new quarks, the 6 6-jet and 26-|-4t final states have already been strongly constrained by the 
current multijet or multilepton searches. Thus in the vector-like quark model, the parameter space for 
dark matter explanation of GeV excess is further constrained. On the other hand, pair production of 
new leptons with 6t final state (with two triple-r resonances, and also resonances in one of the di-r 
pairs in each triplet) is only weakly constrained due to the small production cross section. It will be 
interesting to explore these types of simplified models and new vector-like fermion signatures further 
in the upcoming LHC runs. 
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